Effects of Rashba and Dresselhaus spin-orbit couplings 
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We investigate the pair condensation of a two-spin-component Fermi gas in the presence of both 
Rashba and Dresselhaus spin-orbit couplings. We calculate the condensate fraction in the BCS- 
BEC crossover both in two and in three dimensions by taking into account singlet and triplet 
pairings. These quantities are studied by varying the spin-orbit interaction from the case with the 
only Rashba to the equal-Rashba-Dresselhaus one. We find that, by mixing the two couplings, the 
singlet pairing decreases while the triplet pairing is suppressed in the BCS regime and increased 
in the BEC regime, both in two and three dimensions. At fixed spin-orbital strength, the greatest 
total condensate fraction is obtained when only one coupling (only Rashba or only Dresselhaus) is 
present. 
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I. INTRODUCTION 

In the last years the predicted crossover [H-IH from the 
Bardccn-Cooper-Schricffcr (BCS) state of weakly bound 
Fermi pairs to the Bose-Einstein condensate (BEC) of 
molecular dimers has been observed by several exper- 
imental groups 043) ■ in particular, two seminal ex- 
periments 0, [3> with two hyperfine component Fermi 
vapours of 6 Li atoms in the BCS-BEC crossover, have 
been performed to study the condensate fraction of 
Cooper pairs [10| . which is directly related to the off- 
diagonal-long-range order of the two-body density ma- 
trix of fermions [ill . fl2| . At very low temperature these 
experimental results show an excellent agreement with 
the zero-temperature theoretical predictions of mean- 
field approaches [H, [TiJ and Monte-Carlo simulations 
[15( 1 . However, as discussed in when the effects of 
temperature cannot be neglected it is necessary to in- 
clude beyond mean-field corrections to reproduce quan- 
titatively the experimental data. Intensive theoretical 
studies have been developed on the condensate fraction 
along BCS-BEC crossover for a two-dimensional (2D) 
Fermi gas [l7l - l22l | , and for a three-spin-component Fermi 
gas with SU(3) symmetry f23T - |25j as well. The recent ex- 
perimental realization of 2D degenerate Fermi gases for 
ultra-cold atoms in a highly anisotropic disk-shaped po- 
tential [26] is one of the reasons of the growing interest 
for fermions in reduced dimensionality. 

Very recently artificial spin-orbit coupling has been re- 
alized in neutral bosonic systems [27j . In such systems 
the strength of the coupling can be optically tuned and 
this is indeed a useful tool also for ultracold fermions 
[28l [29j . These achievements have stimulated theoret- 
ical efforts in understanding the spin-orbit effects with 
Rashba [30jj and Dresselhaus [3l| terms in the BCS- 
BEC crossover [H-E3. The evolution from BCS to 
BEC superfluidity was intensively studied in the pres- 
ence of spin-orbit coupling for a 3D uniform Fermi ga s 
ilffl p. El, El and in the 2D case [H, El El® E3 . 
In Ref. [30| we have analyzed the effects of the spin-orbit 



coupling on the condensate fraction by studying both the 
singlet and the triplet pairing contributions in the pres- 
ence of the Rashba coupling. 

In this paper we extend and complete the zero- 
temperature study presented in [39]. In particular we 
investigate the condensate fraction along the BCS-BEC 
crossover with Rashba and Dresselhaus spin-orbit cou- 
plings both in 2D and in 3D analyzing singlet and triplet 
contributions to pairing condensation. We study these 
quantities at zero temperature by gradually including 
the effect of the Dresselhaus term in the spin-orbit cou- 
pling. We show that along this tuning the singlet con- 
tribution to the condensate fraction always decreases, 
while the triplet contribution is strongly suppressed in 
the BCS regime and enhanced in the BEC one, both in 
2D and in 3D. Indeed this enhancement in 2D takes place 
when the binding energy is greater than the Fermi energy 
and in 3D when the dimensionless interaction parameter 
y = l/(kpa s ) (fci? is the Fermi linear momentum and a s 
the interatomic s-wave scattering length) is mainly posi- 
tive. We find that the total condensate fraction is greater 
when only the Rashba coupling is active. Instead, when 
the Rashba and Dresselhaus couplings are equal the total 
condensate fraction is the same as that obtained in the 
absence of spin-orbit coupling. The chemical potential 
and the pairing gap decrease in both the two regimes, 
when the coupling is changed from the only-Rashba to 
equal-Rashba-Dresselhaus case. 



II. THE MODEL 

We describe a gas of tw o-sp in-component Fermi atoms 
with spin-orbit couplings [30L |3l| by using the following 
one-body Hamiltonian 

rfr 2 fc 2 

Hq = ^J^HkVj-^— - + h[v R (a x ky - a y k x ) 

k 

+ V£)(<T x ky + 0~yk X )]W(k), (1) 
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where \x is the chemical potential and vd and vr are, re- 
spectively, the Rashba and Dresselhaus velocities; a x and 
a y denote the Pauli matrices in the x and y directions, 
and ip(k) — (ip^(k),ipi(k)) T is the Nambu spinor. Notice 
that we can rewrite Eq. (Q} as 



mv 



^(k) 



where 



ud 



A = —— 




(sin 9 — cos 6)a y 
sin# + cos 9)a x 




(2) 

(3) 
(4) 

(5) 



In addition to the one-body Hamiltonian Hq we consider 
the two-body interaction Hamiltonian given by 

H! = -i- J2 ^(k + q)^(-k)^(-k' + q)^ t (k') , (6) 

kk'q 

where g > 0, which corresponds to attractive interaction. 
The total Hamiltonian thus reads 



H = H +Hj . 



(7) 



From this Hamiltonian we calculate the effects of Rashba 
and Dresselhaus spin-orbit couplings on singlet and 
triplet condensation with Fermi atoms in the full BCS- 
BEC crossover. 



with 



and 



7(k) = hv R (k y + ik x ) + hv D (k y - ik x ) (14) 



£ k = h 2 k 2 /2m - fx . 
Moreover, the number of particles reads 

£k-|7(k)| £k + l7(k) 



(15) 
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2E 1 (k) 2E 2 (k) 



} , (16) 



while the energy gap A is obtained by solving the corre- 
sponding gap equation 



- = -y(— 



1 



£i(k) E 2 (k) 



(17) 



In addition, the condensate number Nq [52j of Cooper 
pairs is given by 



where 



N C = N + N U 



AT = ^|^ t (k)^(-k))| : 



(18) 



(19) 



1 a L p,f 



16 ^V^i(k) E 2 (k) 



III. GAP ORDER PARAMETER AND 
CONDENSATES 

By decoupling at the mean-field level the two-body 
interaction © we get 

Hl = V l -^f-J2 (A>;(-k)V t (k) + AV|(k)v](-k)) , 

- 1 k 

(8) 

where V is the volume and 

A = (g/V)Y,(M~VMV) (9) 

k 

is the familiar gap order parameter describing the corre- 
lation energy of singlet Cooper pairs. 

As discussed in [39j , from the Hamiltonian ([7]) with Ho 
given by Eq. ((2J and Hi by Eq. ([8]) we can calculate the 
spectrum of single-particle elementary excitations, which 
is given by 

E 1 (k) = y /(e k -|7(k)|) 2 +W (10) 
E 2 (k) = ^ k + | 7 ( k )|)2 + |Ap (11) 
E 3 (k) = -^(k) (12) 
E A (k) = -E 2 (k) (13) 



is the singlet, with total spin 0, contribution to the con- 
densate, whereas 

iVi^K^k^-k)}! 2 (20) 

k 

= w!vf^ i_y 

16 ^V^i(k) E 2 (k)J 

is the triplet one, with total spin 1. Eqs. ([16]), ([17]), ([19]) 
and (|20[) are the starting point of our present investiga- 
tion. The finite-temperature version of these equations 
can be found in our previous paper (39j | . It is important 
to notice that, even if at the mean-field level only the 
singlet energy gap A = Aj^ appears, while the triplet 
one, A<|~|~ = (g/V) 5TJk('/'t( — k)^(k)), is absent, as one 
can see from Eq. ([B]), triplet pairing can be generated by 
the presence of the spin-orbital interaction. 

We study our system both in two dimensional (2D) 
and in three dimensional (3D) case in the absence of the 
temperature, T ~ 0. With this purpose, we have to 
analyze the key quantities for the 2D case - the binding 
energy cb, the condensate densities ns {S = for the 
singlet and S = 1 for the triplet), the chemical potential 
H, and the gap A - and for the 3D case - the interaction 
parameter y, the condensate densities ns, the chemical 
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potential /i, and the gap A. We follow the same path as 
in (39j . that is we define the dimensionless parameters 



x ° = A 



x\ = 2m 
x 2 = 2m 



{vr - vp) 2 _ 2mv 2 
A 

(VR + V D ) 2 



A 

2mv 2 
A A 

so that for the 2D case we have 

2A 2 



(_B_ 

A 



lim 



8 J^-(a;o,xi,X2) 
2ttxq 



a^oo exp[I g (x ,Xi,x 2 )/-K] - 1 
2ns _ 1 i n s (xo,xi,x 2 ) 
n 

f_ = 

(F lff(x ,Xi,X 2 ) 

A 2tt 
e F 

(vr Tv d ) 2 



I 2 j(x ,X 1 ,X 2 ) 



Xl.2 



2 I%(xq,Xx,x 2 ) 



(21) 

(1- sin 26) (22) 
(1 + sin 26>) (23) 



(24) 

(25) 
(26) 
(27) 
(28) 



with A the ultraviolet momentum cut-off. For the 3D, 
one has 



Ia B (X ,X1,X 2 ) 



y ~k F a s 3V37T5/3 J*(a; , a*, a; 2 )V3 

2n s 1 Iff s (xo,xi, x 2 ) 
n 8 I%{x ,xi,x 2 ) 



/ nr \ 2/3 

a;o^Ar(a;o,a;i,a;2) _2/3 



ii 

A 



= 4 



(|) /3 /at(o;o,xi,X2) 2/3 



(29) 
(30) 
(31) 
(32) 



(^) V3 xi, 2 I%(xo,x 1 ,x 2 )- 2/3 (33) 



with fci? the Fermi wave vector and a s the interatomic 
s-wave scattering length. All the above equations are 
written in terms of the following dimensionless integrals: 



1 f A 

I g (x ,xi,x 2 ) = - / d 2 q ^2 



lff(x ,xi,x 2 ) = / d d q 



1 \j {q 2 - x + r x x ql + x 2 ql)j +1 
^ <? 2 - x + rJxiql + x 2q 2 ^ 



r ±x \J (q 2 - Xo+r^J xi ql + x 2 q^j + 1 y 



E 



' 1 ■ ' ' q 2 — x + r x lx 



iq 2 + x 2 qfj +lj 



1 1 

^2 ~ 2 



E 



r 



(34) 



where d = 2 holds for the 2D case, while d = 3 for the 
3D one; S = for the singlet and S = 1 for the triplet 
contribution. 



IV. RESULTS 

In Fig. [T] we plot the condensate fraction of the 2D 
Fermi system as a function of the scaled binding energy 
£s/e_F, with cf the 2D Fermi energy, for different values 



of the characteristic velocity v 



and two val- 



ues of the mixing angle 9 = wct&n (vd/vr). In the first 
two panels there are the singlet (top panel) and triplet 
(middle panel) contributions to the condensate fraction, 
while in the lower panel there is the total condensate 
fraction. The results are shown for 9 = (vd = 0, solid 
curves) and 9 = tt/4 (vr = Vd, dashed curves) and three 
values of the velocity v. 

In Fig. [2]we plot instead the condensate fraction of the 
3D Fermi system as a function of the scaled interaction 
strength y = l/(kpa s ) with kp the 3D Fermi linear mo- 
mentum and a s the s-wave scattering length. As in Fig. 
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FIG. 1: (Color online) 2D Fermi superfluid. Singlet (top 
panel) and triplet (middle panel) contributions to the con- 
densate fraction for v 2 = 1, 2, 4 (in units of v F ) and for 6 = 0, 
i.e. only Rashba term, (red solid curves) and 9 = 7r/4, i.e. 
for equal Rashba and Dresselhaus terms, (blue dashed line). 
Bottom panel: the total condensate fraction, i.e. the sum of 
both singlet and triplet contributions, for the same values of 
v and 9. Notice that for 8 — 7r/4 (blue dashed line), the full 
condensate fraction is the same as without spin-orbit (v = 0). 



FIG. 2: (Color online) 3D Fermi superfluid. Singlet (top 
panel) and triplet (middle panel) contributions to the con- 
densate fraction for v 2 = 1, 2, 4 (in units of Vp) and for 8 — 0, 
i.e. only Rashba term, (red solid curves) and 8 — 7r/4, i.e. 
for equal Rashba and Dresselhaus terms, (blue dashed line). 
Bottom panel: the total condensate fraction, i.e. the sum of 
both singlet and triplet contributions, for the same values of 
v and 9. Notice that for 9 = 7r/4 (blue dashed line), the full 
condensate fraction is the same as without spin-orbit (v = 0). 



Q] in the first two panels there are the singlet (top panel) 
and triplet (middle panel) contributions to the conden- 
sate fraction, while in the lower panel there is the total 
condensate fraction. Again, results are shown for 9 = 
(vd = 0, solid curves) and 6 — 7r/4 (vr = Vn, dashed 
curves) and three values of the velocity v. The two sets 
of plots corresponding to 2D and 3D share many features 
and would have looked very similar if we had plotted Fig. 
Q]in terms of ln(l/fci?a2D) = ln(e.B/2eF)/2. 

Our calculations show that the singlet contribution 
to the condensate fraction, 2no/n, both in 2D and in 
3D, decreases when one fixes v/vp and moves from the 
only-Rashba (or only-Dresselhaus) case (9 = 0) to equal- 
Rashba-Dresselhaus case (8 = ir/4), as it can be ob- 
served from the top panels of Figs. [Hand [2] In the case 
of only Rashba (9 — 0) we observed [39[ that, for one 



value of the scattering parameter close to the crossover 
(tB — £f in 2D and y ~ in 3D) and for v > vp, 
no/n does not depends on v, namely there is a nodal 
point for the singlet condensate fraction when one in- 
creases largely enough the Rashba coupling, see red solid 
curves in the top panels of Figs. [T] and [21 Regarding 
the triplet contribution, 2ni/n, we find that for suffi- 
ciently small values of the binding energy £b/^f (in 2D) 
or for y < (in 3D), 2ni/n is suppressed by mixing 
the two spin-orbital couplings, whereas when cb/^f (in 
2D) becomes large enough (eb > 1.5ejr) and y > (in 
3D), the triplet contribution is enhanced, see the middle 
panels of Figs. [TJ and [2j The full condensate fraction, 
nc/n = 2(uq + n{)/n, at fixed v/vf, is maximum for 
only-Rashba (or only-Dresselhaus) case (9 = 0). In the 
extreme case of equal Rashba and Dresselhaus contribu- 
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/i ~ €f — mv 2 . In the case of 9 = 7r/4, Eq. (j36|) is valid 
for all values of the binding energy cb and comes simply 
from gauge transforming the fields, ip(r) — > e lyA «^p(r), in 
the Hamiltonian Eq. @. 

Let us now investigate in detail the condensate frac- 
tions no/n and ni/n, the chemical potential /i and the 
energy gap A when 9 is increased from to 7r/4. For com- 
putational simplicity, we consider the two dimensional 
case. 
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FIG. 3: (Color online) 2D Fermi superfluid. Top panel: Sin- 
glet (red solid line, left y-axis) and triplet (green dashed 
line, right y-axis) contributions to the condensate fraction, 
for v 2 = 2vp and es = 0.5eF, as functions of 9. Notice that 
2no/n and 2ni/n are plotted with different scales. Middle 
panel: chemical potential, fj, (in units of ep), for v 2 — 2vp 
and 6b = 0.5eF, as a function of 9. Bottom panel: gap func- 
tion, A (in units of ef), for v 2 = 2vp and es = 0.5eF, as a 
function of 6. 




tions (9 = 7r/4), instead, the total condensate fraction is 
the same as that obtained without spin-orbit at all, i.e. 
with v = 0, see the dotted line of the bottom panels of 
Figs. [T]and[2] In particular, for the case of 9 = 7r/4 the 
following results for the condensate fraction, the chemical 
potential and the gap function hold 

— M = tt/4) = —(v = 0) (35) 
n n 

H(v,9 = ir/4) = fi(v = 0)-mv 2 (36) 

A(u, 6 = tt/4) = A(v = 0) . (37) 

The chemical potential is consistent with the known 
perturbative result for small v and small es, which is 



FIG. 4: (Color online) 2D Fermi superfluid. Top panel: Sin- 
glet (red solid line, left y-axis) and triplet (green dashed 
line, right y-axis) contributions to the condensate fraction, 
for v 2 — 2vp and es = 2ef, as functions of 9. Middle 
panel: chemical potential /i (in units of ep) for v 2 — 2vp 
and es = 2ep, as a function of 9. Bottom panel: gap func- 
tion, A (in units of ef), for v 2 = 2vp and es = 2ep, as a 
function of 9. 

In Figs. [3] and S] we plot these quantities by choosing 
v 2 = 2vp in correspondence to two different values of 
the binding energy. These two binding energies are those 
for which a softening and a hardening of the triplet con- 
tribution (2m/n) against 9 is expected; thus Fig. [3] is 
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obtained with cb = 0.5£f and Fig. [4] with e B = 2ep (see 
the above discussion about Fig. [IJ. From the former fig- 
ure, Fig. it can be pointed out that, for small binding 
energies, Inyjn decreases by increasing 9 from to 7r/4, 
and the same behavior, but rescaled, is observed for the 
singlet contribution 2no/n (see the top panel). Both the 
chemical potential and the gap function decrease going 
from the only-Rashba case to a fully mixed one. Fig. SI 
instead, shows that, for large binding energies, 2n\/n in- 
creases when 9 is changed from 9 = to 9 — tt / 4, while 
2no/n always decreases, as well as the chemical potential 
and the gap function. In both above discussed cases, the 
singlet contribution to the condensate fraction is always 
greater than the triplet one. 

V. CONCLUSIONS 

We have analyzed the condensation of fermionic atoms 
along the BCS-BEC crossover in the presence of Rashba 
and Dresselhaus spin-orbit couplings. The condensation 
has been characterized by calculating the singlet and the 



triplet contributions to the pairing and therefore the full 
condensate fraction. We have studied these quantities 
by varying the spin-orbit from the situation in which the 
only Rashba coupling is present to that in which the 
Rashba and Dresselhaus velocities are equal. We have 
found that moving along this path, the singlet contri- 
bution to the condensate fraction decreases, while the 
triplet one behaves differently in the two regimes (BCS 
and BEC). In the BCS regime, the triplet pairing is sup- 
pressed upon mixing the two spin-orbit couplings while 
in the BEC regime it experiences an enhancement over 
the only-Rashba case. In other words, the triplet pairing 
is maximized in the BCS regime if only Rashba (or only 
Dresselhaus) term is active, while it is strengthened in 
the BEC regime by mixing the two spin-orbital couplings. 
This behavior takes place both in two and three dimen- 
sions. We have shown also that the total condensate frac- 
tion is greater when only one coupling (only Rashba or 
only Dresselhaus) is present, while in the equal-Rashba- 
Dresselhaus case, is the same as that obtained without 
spin-orbit. 
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